3D XY scaling theory of the superconducting phase transition 
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The intermediate 3D XY scaling theory of superconductivity at zero and nonzero magnetic fields 
is developed, based only upon the dimensional hypothesis B ~ (Length) -2 . Universal as well as 
nonuniversal aspects of the theory are identified, including background terms and demagnetization 
effects. Two scaling regions are predicted: an "inner" region (very near the zero field supercon- 
ducting transition, T c ), where the fields B, H, and H cyi differ substantially, due to the presence of 
diamagnetic fluctuations, and an "outer" region (away from T c ), where the fields can all be treated 
similarly. The characteristic field (-ffo) and temperature (ti) scales, separating the two regimes, are 
estimated. Scaling theories of the phase transition line, magnetization, specific heat, and conductiv- 
ity are discussed. Multicritical behavior, involving critical glass fluctuations, is investigated along 
the transition line, T m (B), at nonzero fields. 

74.40.+k,74.25.Bt,74.25.Dw; Keywords: critical phenomena, fluctuation effects, mixed state, phase 
diagram. 



I. INTRODUCTION 



The phase transition into the superconducting state in 
the high-T c oxide materials has been the focus of much 
theoretical and experimental research. This transition is 
particularly interesting because mean field theory does 
not provide an adequate description; for both zero and 
nonzero magnetic fields, fluctuations play an important 
role. The challenge is to acquire a sound understanding 
of these fluctuations in different regimes of interest. 

The superconducting transition corresponds to the 
melting of a vortex solid, and occurs along the line T m (B) 
in the field-temperature plane. When B > 0, the tran- 
sition can be either first order, as seems to be the case 
for pure systems or continuous, when systems con- 

tain disorder Q. However, we must distinguish between 
the transitions at zero and nonzero fields. For example, 
although T m (B > 0) may be first order, the zero field 
transition, T c , is expected to be continuous. The zero and 
finite field transitions must therefore belong to different 
universality classes, as corroborated experimentally [HQ. 
If the transition T m (B > 0) is also continuous, it must 
be associated with a second type of critical behavior. At 
low fields, it is possible to observe both critical behaviors 
simultaneously ■ The fact that the B > transition 
line, with its distinct critical behavior, connects smoothly 
to the zero field transition [T c = lime— »o T m (B)] implies 
that (T = T C ,B = 0) is a multicritical point ||. 

As happens in two dimensional superconductors, the 
three dimensional (3D) superconducting phase transition 
occurs somewhat below the temperature where fluctua- 
tions of the order parameter amplitude drive the density 
of fluctuating Cooper pairs to zero. A growing body of 
experimental evidence [|]]|,||-||,|l(^-|l7| supports the idea 
that the zero field transition of strongly type-II supercon- 
ductors, such as the high-T c oxides, is in the "intermedi- 
ate" 3D XY model class, which includes the A-transition 



in 4 He. The essential fluctuations of this class involve the 
order parameter phase, but do not involve its amplitude 
or the vector potential. 

When vector potential fluctuations are also considered, 
in addition to phase fluctuations of the order parame- 
ter, the ultimate critical behavior of the superconducting 
phase transition may belong to the "inverted" 3D XY 
model class The relative stiffness of the vector po- 
tential in strongly type-II superconductors, compared to 
the order parameter phase, causes this inverted critical 
regime to be very small. In this work, we focus our at- 
tention on the intermediate scaling region. 

Critical fluctuations along the transition line T m (B > 
0) are thought to be of the glass type, in the case of strong 
disorder Q,|l9|,^0|] . However, experiments show that 
3D XY fluctuations continue to be relevant for the scal- 
ing when B > 0. These include measurements of the spe- 
cific heat H| 10 OJl^l, magnetization [Bp II -U|, pen- 
etration depth |1J, and conductivity ||"i'|lj,|17[. While 
3D XY and glass fluctuations may coexist at low fields 
near the superconducting transition, they correspond to 
two distinct universality classes, with distinct exponents 
and scaling functions. This multicritical coexistence of 
superconducting fluctuations has been studied experi- 
mentally ||-|. 

In this paper we consider intermediate 3D XY fluctu- 
ations, for low magnetic fields, B > pi}] . For simplic- 
ity we consider only the case of a continuous transition 
T m (B) at nonzero fields; details of the case of first order 
melting are presented elsewhere 22 . Some other theoret- 
ical treatments of the 3D XY model of superconductivity 
include Refs. [ll^Ji"9| , p3| -p(| , while some relevant numer- 
ical simulations are given in Ref. |27]]. Although exact 
solutions to this problem are not yet in sight, a fruit- 
ful advance can be made using the scaling approach. In 
this method, a scaling form for the free energy is hypoth- 
esized by means of a dimensional analysis. Using this 
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ansatz, the theory is developed in a very general form, 
relying as little as possible on particular models for the 
superconductor. In addition to the free energy, we dis- 
cuss the magnetization, specific heat, and conductivity. 
One important aspect of the present work is to clearly 
specify the nonunivcrsal parameters which enter the the- 
ory. These are of three types. One type comes from the 
smooth background terms. The second type is associated 
with demagnetization effects, and is related to sample 
geometry. The third type involves material-dependent 
constants which enter the scaling term of the free energy. 

We also carefully distinguish between the different 
fields: the spatially averaged magnetic field B, its con- 
jugate field H, and the external field H cyi . The usual 
assumption that B ~ iJ ex , appropriate for a disk geome- 
try, is reconsidered. The distinct thermodynamic roles of 
B and H, defined through the relation H = 4ndf/dB (f 
is the appropriate free energy density), suggest that the 
two quantities should scale differently; this reflects the 
emergence of diamagnetic fluctuations. The distinction 
between the different fields is most apparent in a small 
"inner" scaling region near the zero field transition. In 
a realistic physical scenario, H cx is the externally con- 
trolled variable; B and H then both acquire fluctuation 
contributions. Since estimates of the inverted XY scaling 
regime place it in the vicinity of the inner scaling region, 
it is therefore crucial to elucidate the differences between 
the different fields, in order to unravel the different crit- 
ical phenomena. 

The plan of the paper is as follows. Section II dis- 
cusses the scaling form of the free energy, Eq. (^). Sec- 
tion III derives the magnetic equation of state, Eq. (||). 
The crossover field, Hq, between the inner and outer scal- 
ing regions, is identified in Eq. (||). Section IV derives the 
form of the superconducting phase transition line in the 
B-T, H-T, and H cx -T planes, given by Eqs. (§), 
and ( |l4|) , respectively. The H-T phase diagram is shown 
in Fig. 1. Section V is a brief interlude, which shows how 
the Abrikosov theory of the superconducting transition 
in a field is a special case of the more general scaling the- 
ory. The magnetization is discussed in Section VI, where 
we derive the relations between M and H [Eq. (^5|)] or 
M and H cx [Eq. (||)], for T = T c or T m . In Section VII 
we obtain the specific heat, Eq. (p9|). Section VIII pos- 
tulates the dynamic scaling theory associated with the 
ohmic conductivity, Eq. (|32]). In Section IX we conclude, 
giving estimates for the size of the inner and inverted 
scaling regions. 

Notation: In this paper, only the quantities with tildes 
involve factors related to sample geometry. 

II. SCALING THEORY 

We now discuss the basic thermodynamics of the scal- 
ing approach. At zero field, we adopt the usual scaling 



hypothesis, which states that any observed singular be- 
havior involves the divergence of the correlation length, 
in terms of the relative temperature t = (T—T c ) /T c . The 
zero field correlation length £(T) and specific heat C(T) 
are then described by power laws, with the exponents v xy 
and a xy , respectively pS} ]. 

To extend scaling to finite fields, we must know how 
B scales. We therefore introduce a definite assumption 
about the physics of a superconductor. A defining char- 
acteristic of a superconductor is its broken U(l) or gauge 
symmetry, which is reflected in the symmetry of the su- 
perconducting order parameter. Gauge-invariance then 
implies the following identification for the gradient op- 
erator: V — ► D = V + 2ieA/hc. The basic scaling ar- 
gument, which amounts to a dimensional analysis, states 
that the two terms appearing in D must have the same 
scaling dimension: (Length) . Similar dimensional ar- 
guments were proposed in Ref. Jl9[ , and later confirmed 
in Ref. p6| . The dimensionality of the magnetic field is 
then expressed as 

B = |V X A| ~ (Length)" 2 . (1) 

The most general scaling hypothesis for the free energy 
density becomes Q 

f{B,T) = f b {B,T) + f k \t\ 2 - a * y <t>± ( ^hI^ )- ( 2 ) 

The function fb(B, T) represents the smooth background, 
while the second term in Eq. (Q) encapsulates all the ef- 
fects of 3D XY critical fluctuations. The scaling theory is 
universal in the sense that neither the exponents, v xy and 
a xy , nor the functions 4>+ {4>-)i corresponding to t > 
(t < 0), contain any sample dependence; they are the 
same for all superconductors which exhibit 3D XY scal- 
ing. The sample dependence rests only in the background 
term fb{B, T), the transition temperature T c , and in the 
parameters fk and Hk in the fluctuation term. Here, fk 
has units of free energy density, making <p± (x) dimension- 
less, while Hk has field units, making the scaling variable, 
x = B\t\~ 2Vscy / Hk, dimensionless. 

We now discuss several points concerning the free en- 
ergy, Eq. (|): 

(1) The arguments leading to Eq. (Q) do not admit 
anomalous scaling dimensions, since the magnetic field 
does not fluctuate in the intermediate XY model class. 
On the other hand, fluctuations of the vector potential 
become important in the inverted XY critical region, 
and may lead to the appearance of an anomalous di- 
mension: B ~ (Length) ~ 2+1? . Mounting experimental 
evidence seems to be consistent with the intermediate 
scaling presciption of § = Q. In the work which follows, 
the magnetic field is always treated within the intermedi- 
ate scaling hypothesis. In particular, we emphasize that 
the inner scaling region, discussed below, also emerges 
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from Eq. (0), and is not related to inverted 3D XY be- 
havior. 

(2) The region of the B-T plane for which Eq. (||) pro- 
vides the correct scaling description must either be de- 
termined experimentally, or by a more detailed theory. 
Empirically, the range of validity of Eq. (||) appears to 
extend far beyond the characteristic 3D XY temperature 
and field scales, t\ and Ho (discussed below), which are 
very small in many materials. 

(3) Near T — T c , the background term in / may be 
written approximately as follows: 

f b (B,T)~f (T) + f 2 (T)B 2 , (3) 

where fo(T) and f 2 (T) are smooth functions, with no 
singular behavior. We shall see that f 2 {T) is related to 
the background magnetic susceptibility. 

(4) The sample dependent quantities defined above, 
T c , fk, Hk, and f b (B,T), must all reflect the anistropy 
of the superconductor. The main effect of anisotropy on 
Eq. (||) is that H^ 1 becomes a tensor quantity. In other 
words, anisotropy causes Hk to depend on the direction 
of B. To avoid such complications here, we may simplify 
the analysis by choosing our geometry carefully: we con- 
sider the simple, but experimentally relevant case that 
anisotropy, if present, is aligned with the principal axes 
of an assumed ellipsoidal sample. In addition, we require 
that any external magnetic field should be applied along 
a principal axis of the sample. In other geometries, the 
field B becomes nonuniform and/or misaligned with the 
external field H ox . 

(5) It is possible to relate the parameter fk to an- 
other nonuniversal parameter £oj which appears in the 
zero field critical correlation length £(T) = Col^l - ^" , by 
means of two-scale- factor universality |^5],|30| . The result 
can be written as fk = /cbT c /£q, where for simplicity, we 
have absorbed a universal proportionality constant into 
the definitions of fk and £o- Without loss of generality, 
we may also introduce anisotropy into this relation as 
fk = jk B T c /£l w , where 7 = £ a &o/£cO represents the ra- 
tio of the zero field 3D XY critical correlation lengths 
along different axes. 

Several theoretical analyses of superconductors re- 
late the renormalized anisotropy parameter, 7, to the 
anisotropy appearing in the bare Hamiltonian, by show- 
ing that the anisotropic problem may be treated as an 
isotropic one |]l6| , ^l| . In this paper, we also assume 
that the isotropic and anisotropic problems should in- 
volve the same, universal scaling functions. However, we 
work strictly in terms of the measurable (renormalized) 
anisotropy factor. Note that unless ^ a &o can be deter- 
mined by independent means, it is not possible to extract 
7 directly from Eq. (^J), without making further assump- 
tions. (These are described in the following point.) We 
therefore retain the more general (fk, Hk) notation here, 
noting that anisotropy is naturally absorbed into these 
parameters. 



(6) The scaling ansatz of Eq. (||) is very general, since it 
arises from purely dimensional arguments. This form in- 
volves exactly two nonuniversal parameters ]3C| ], fk and 
Hk, in addition to the temperature scale, T c , and the 
background term, f b (B,T). It is not possible to reduce 
this number of sample dependent parameters without 
further information. Furthermore, we emphasize that the 
3D XY model has not been solved exactly, and there- 
fore cannot provide such information. However, it has 
been suggested that a relation exists between fk and Hk 
p9| , p5| , thus reducing the number of nonuniversal param- 
eters by one. 

The heuristic argument states that a characteristic 
field scale, B c h(T), appears in the argument of the scal- 
ing functions 4>±(x), in the form x = B/B c h(T). This 
field scale should be given precisely by [jl9|,|5| B c h(T) — 
§o/£, 2 {T), where £(T) is again the zero field correla- 
tion length. It follows that Hk — ^o/foi fr° m which 
we obtain the desired relation: /fc/fcsT c = (Hk/^o) 3 ^ 2 ■ 
For anisotropic superconductors, the relation becomes 
fk/jk B T c = (iJ fe / $ o) 3/2 , when B||c, thereby providing 
a method for determining the anisotropy, 7, when the 
parameters fk and Hk are determined experimentally. 

The proposed relation between fk and Hk can be tested 
through a scaling analysis of the fluctuation mangetiza- 
tion, as described in Section VI. We point out that the 
relation between certainly does not hold for the case of 
the mean field Abrikosov theory, described in Section V. 
However, for critical fluctuations, the question is still 
open. For generality here, we continue to treat fk and 
Hk as independent parameters. 

III. EQUATION OF STATE 

To procede with the analysis of Eq. ^ , the hyperscal- 
ing relation may be used: 2 — a xy = Zv xy (for the 3D 
case). Eqs. (||) and (^) can be rewritten as 

f(B,T) = f (T) + f 2 (T)B 2 + fk\t\ 3 ^^ ± ( m H 2 J* V ^ ■ 

(4) 

The magnetic equation of state is derived from the iden- 
tity H = Andf/dB: 

(5, 

where (f>'± (x) = d(fi±/dx. Eq. (|5|) relates B and H for 
a given temperature. The normal and fluctuation con- 
tributions appear in the first and second terms, respec- 
tively. Ht = 871-/2 (T) is the inverse permeability of the 
background; experimentally, it is found that Qt — 1- 

It is desirable to invert Eq. (||) to obtain B(H,T). In 
general, this is not possible, because 4>±{x) are not yet 
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known theoretically. However, progress can be made in 
certain cases. When B > and T = T c , the fluctuation 
part of Eq. (0) must be smooth and independent of t. 
This leads to the following asymptotic behavior: 



lim 4ir(f>' ± (x) = (b x) 1/2 , 



(6) 



where, bo is a dimensionless, universal constant of the 
3D XY theory. Thus when T = T c , we have 



n c B = H + 2H Q - 2H„ /2 V / H +H. 

Here we have simplified the notation using il c = Qt c . 
In Eq. (j^) we have defined a characteristic field 

H = bofl/4ft c Hl 



(7) 



(8) 



which represents the crossover between the inner and 
outer scaling regions, for T — T c . This crossover field 
is very small for most superconductors, due to the small 
size of diamagnetic fluctuations above the transition. In 
the outer scaling region, H ^> H , we see that the usual 
approximation H ~ Q C B is appropriate. However, in 
the inner region, H <C Ho, diamagnetic fluctuations 
dominate the equation of state, leading to the behavior 
B oc H 1 / 2 , for T = T c . If there exists an exact rela- 
tion between and Hk, as described in Sec. II, then 
Eq. (||) should reduce to H oc j 2 T 2 , demonstrating the 
strong anisotropy dependence of the inner scaling region 
becomes explicit. 



nearly constant over the temperature range of interest: 
Eq. ( |l0| ) then becomes 



H(T) 



(11) 



^cXm^-k 

We have introduced the relative temperature scale 

h = (/^AW^ 2 ) 1 ^, (12) 

and the universal number b\ = 4 7r<//_ (£,„). The expres- 
sion t\T c represents the size of the inner scaling region 
along the transition line T m (H); it is analagous to Ho 
on the H a> ase diagi in 

Fig. 1. 
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IV. PHASE TRANSITION LINE 



The free energy (^) must contain information concern- 
ing the phase transition into the superconducting state. 
Of crucial importance is the fact that the superconduct- 
ing transition line T m (B) terminates at T c on the B = 
axis. The special point (T = T c , B — 0) is then mul- 
ticritical For the case of a continuous transition 

considered in this paper, the transition always occurs 
at a particular, universal value of the scaling variable: 
x = x m . The transition line is then given by 



FIG. 1. Phase diagram in the H-T plane, showing 
the phase transition line T m (H) and the different scal- 
ing regions. In the shaded, "inner" scaling region, the 
usual approximation B ~ H breaks down. The dimen- 
sions of the inner region are of order Ho x tiT c , where 
H is defined in Eq. (|) and ti in Eq. ([jj). In the 
unshaded, "outer" region, the approximation B ~ H 
is accurate. Relative sizes of the two scaling regions 
are not drawn to scale; in many cases, the inner region 
is tiny compared to the outer region. 



B(T) - {x m H k )\t\ 



(9) 



Combining Eqs. (H) and (JoJ) , we obtain the transition 
line in the H-T plane: 

H(T) = n T x m H k \t\ 2 ^ + ^J±^( Xm )\t\"*y . (10) 



Eq. (10) has two terms, one going as lij^", and the other 
going as \t\ 2 " my . Two different types of nonuniversal pa- 
rameters enter the \t\ 2V:cv term, including the background 
quantity fix, which contains a possible temperature de- 
pendence. For simplicity, we will assume that fir is 



It is also of interest to study the phase transition line in 
the H cx -T plane, because of its experimental importance. 
We assume that H cx is applied parallel to a principal axis 
of our sample, which is an ellipsoid [ 33 1 . A simple equa- 
tion then relates the different fields [ 34 1: 



H ex =nB + (l- n)H, 



(13) 



where n is the demagnetizing coefficient, satisfying < 
n < 1. (n ~ 0.8 — 0.9 for typical high-T c single crystals, 
while n > 0.99 for thin films.) Eqs. ®, ©, and 
then give the transition line: 



4 



(h\t\r 



(14) 



where _(1 - T ) = (1 - n)(l - Q T ), and t" mv = < xy (l - 
n)il c /Ct c . Note again that quantities with tildes differ 
from those without tildes only by geometric factors. 

The limiting behaviors of Eqs. (11) and ([TJ) are clear; 
we show only the result for H cx : 



(h c x m H k t^)\t\^y |t|«t~i, 

(n cXm H k )\t\ 2 ^ |t|»t x . 



(15) 



The large \t\ limit corresponds to the usual (outer region) 
description of the phase boundary ]l9[ . New behavior is 
observed in the small inner region, \t\ <C t±. 



V. MEAN FIELD THEORY 

As noted in the introduction, the superconducting 
transition in high-T c materials should not be amenable to 
mean field description. However in low-T c superconduc- 
tors, the mean field description is often found to be accu- 
rate. It is therefore interesting to note that the Abrikosov 
theory |35) of the B > mean field transition is consis- 
tent with the scaling theory presented so far. We now 
pause briefly to discuss this point. 

Beginning with Eq. (Q), progress cannot be made us- 
ing hyperscaling, since this relation does not hold for the 
mean field transition. However, mean field exponents 
may be used, since they are known exactly: v m f = 1/2 
and a m [ = 0. This gives 



/mf = fo(T) + f" 

8n 



/fc|to| 2 <^mf± 



B\to 



(16) 



where to = (T — T c q)/T c q, and T c q refers to the temper- 
ature where the upper critical field vanishes: B c2 {T = 
T c o) = 0. For simplicity, we have assumed that the nor- 
mal state background has no magnetic effects (fir = 1), 
as is usual in the Abrikosov theory. The equation of state 
then becomes 



H = B 



Hi. 



M0r 



mf± 



B\t \ 



(17) 



The mean field phase transition, B c2 (T), occurs when 
the argument of </> m f_(:r) has the value x m , leading to 
the following transition line in the B-T plane: 



B = (x m H k )\t \. 



(18) 



Using Eq. (|l7|), the transition line in the H-T plane is 
given by 



H 



x m H k + —-—<P m i-{x m ) 



T, 



cO 



T 



T, 



,.-o 



(19) 



The linearity of the transition line in T c q — T is consis- 
tent with the Abrikosov theory. We emphasize that this 
prediction, Eq. (fl9|), has been obtained with no explicit 
knowledge of the Abrikosov (B > 0) solution. Instead, it 
is a general consequence of the scaling theory. 

The scaling functions <p m f±(x) can now be explicitly 
computed. Recall the Abrikosov solution for the free en- 
ergy, f m {(B,T), near the upper critical field H c2 (T) fl35| : 

(H c2 -B) ~ 

/mf , , 

B > H c2 (T) 



3tt 1+(2k 2 -1)/3a 

Ml 

8tt 



where 



H c2 (T) = 



T c0 -T 



27re 2 



BCS 



T, 



cO 



(20) 



(21) 



Here, k is the Ginzburg parameter, Pa — 1-16 for the 
triangular vortex lattice, and £bcs ~ I^vf /ksTcQ is the 
temperature independent coherence length. 

The correspondence between Eqs. ( |l6| ) and ( p0| ) be- 
comes transparent by making the following identifica- 
tions. The nonuniversal parameters can be taken as 
H k = $o/2<l cs and f k = H 2 k [l + (2 K 2 + l)^}' 1 /Sir, 
where we have adopted the following normalization: 
<?W-(0) = —1. The quantities f k and H k contain all 
the sample dependence of the scaling description. Com- 
parison with Eq. (Qq) now gives 



(pmi-(x) - 
4>mf+(x) = 0. 



-(l-x) 2 




x < 1 
x > 1 



(22) 
(23) 



Within the mean field approach, we can evaluate the uni- 
versal quantities b\ and x m introduced earlier. We find 
that b\ — 0, reflecting the lack of fluctuations in this 
case. (At temperatures above the transition, we simply 
have B = H .) Additionally, we find x m = 1. 

Finally, we point out that mean-field theory provides 
an example of a case where the nonuniversal parameters, 
f k and H k , cannot be related in a simple way, except 
by introducing an additional, nonuniversal parameter k. 
However, H k does take the form suggested by heuristic 
arguments in Sec. II. Note that the relation found here, 
f k oc Hj*, differs from the one described in Sec. II, due to 
the breakdown of hyperscaling in mean field theory. 



VI. MAGNETIZATION 

In the remainder of this paper we consider only 3D XY 
critical fluctuations. As usual, the magnetization is de- 
fined by M = (B- H)/4n. Using Eq. (§|), we obtain 



M(B, T) 



i - n 7 

47T 



-B 



fk 



H 



3/2 



B 1/2 M- 



H~ k 



(24) 



5 



where M.±{x) — x^ 1 / 2 0' ± (a;) are universal scaling func- 
tions. The first term in Eq. ( p4[ ) represents the normal 
state background, and vanishes for Q,t = 1. The second 
term is the critical fluctuation contribution. 

As discussed in Sec. II, certain physical arguments 
may lead to a reduced number of nonuniversal param- 
eters, by providing a relation fkH k 3 ^ 2 cx jT c (7 = 1 for 
the isotropic case), where the proportionality constant is 
not sample dependent. Eq. offers a convenient ex- 
perimental test of this prediction, since the expression 
fkH k , can be directly infcrcd from the scaling. To 
test the prediction, the quantities T c and 7 should be 
determined by independent methods p^] . 

Deriving equations for M(H,T) or M(H ex ,T) is not 
straightforward, because of the difficulty in inverting 
Eq. (|) to obtain B(H, T). We therefore limit our deriva- 
tion to the formulae relating M and H, or M and H ox , 
at the special temperatures T c and T m . Using the appro- 
priate equations of state, we find the following results at 
T = T r : 



4ttO c M = (1 - fl c )H + 2H Q 
4tt(1 - n)Q c M = (1 - Q c )H ex + 2H 



- 2H, 



1/2 



Ho + H c 



(25) 



(26) 



where Hq = Ho(l — n) 2 Vl c /Vl c . The first term on the 
right hand side of both equations gives the normal back- 
ground contribution, while the remaining terms represent 
the fluctuations. The asymptotic behavior of the fluctu- 
ation part of the magnetization, Mfl, can be found. We 
show the results for H PX : 



He: 



M fl 



47r(l— n)f2 e 
27r(l-ri)I2 c 



H cx <g; Hq 

Hex 3> Hq 



(27) 



In the low field limit, Mfl becomes asymptotically linear 
in H ex , like the background. The fluctuation magneti- 
zation then dominates over the background by a factor 
proportional to (I — n) _1 (I — r2 c ) _1 , which can be quite 
large for typical (flat) samples. Thus, in the low field 
(H cx -C Hq), inner scaling region, the problem of back- 
ground subtraction, which otherwise troubles experimen- 
tal analyses, is ameliorated. 

When T = T m , we may still use Eqs. @ and ^ by 
making the following replacements: 



H 



hi 



x m b 



Hr> and Ho 



hi 



x m b 



Hn 



(28) 



VII. SPECIFIC HEAT 



C{B,T) 
T 



= -fo(T) m)B 2 

+ f k T- 2 \t\- a °>y^± 



B\t\ 



Hi- 



(29) 



The dimensionless functions ip±(x) depend on <f)±(x), and 
their first and second derivatives. The first two terms in 
Eq. ( p9|) represent the background, while the last term is 
due to superconducting fluctuations. 
We note the following points: 

(1) Experimentally, it is difficult to isolate the fluctua- 
tion contributions in Eq. (]2£j), due to (i) the smallness of 
fluctuations compared to the background, (ii) the weak- 
ness of the specific heat singularity jj7| (a xy ~ — 0.01), 
(iii) rounding effects, which are often observed in exper- 
iments. 

(2) It is possible to simplify Eq. (H|) for the high-T c 
materials, by using the empirical fact that /^'(T) ~ 0. 
The following scaling quantity may then be considered 



AC(B, T) = C(B, T) - C(0, T), 



(30) 



which contains no background dependence. However, if 
the zero field specific heat cusp becomes rounded near T c , 
as is often the case for real samples, then the imperfect 
scaling, associated with the rounding, is transmitted to 
AC(B,T). 



VIII. CONDUCTIVITY 



Up to this point, we have developed our scaling anal- 
ysis for thermodynamic quantities. All results have been 
derived from the expression for the free energy density, 
Eq. d). However, the description of transport measure- 
ments, such as the conductivity, requires further informa- 
tion. Since the equations of motion governing the time 
evolution of the superconductor are not well understood, 
a dynamic scaling ansatz must be postulated. 

Let us consider the ohmic conductivity a, in order to 
avoid complications arising from current dependence [p8[ . 
For B = 0, when approaching T c from above, a diverges 
according to some power law. Fisher, Fisher, and Huse 
have given arguments leading to the following ansatz [fl9[ : 



<Tfl cx f 



5 = 0, 



(31) 



where cta is the fluctuation part of cr, and z xy is the expo- 
nent of the supposed 3D XY dynamic universality class. 
Scaling at finite fields then proceeds in the usual way: 



The free energy (j|) may be used to compute the spe- 
cific heat at constant B: 



a = S b (B,T) + S k \t\-"°y( z "'- 1) i: ± 



B\t\ 



Hi. 



(32) 
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The first term, Sb{B,T), represents the smooth back- 
ground conductivity. The second term is the fluctuation 
contribution, where the sample-dependent parameter S k 
has dimensions of conductivity. S+(aj) [£_(&)] should 
be universal scaling functions, corresponding to t > 
[t < 0]. The background conductivity term plays a rel- 
atively small role near the transition line, T = T m (B), 
due to the divergence of aa- 

The functions E±(x) are not known theoretically, but 
we may deduce their asymptotic behavior. For B > 
and T — > T c , the conductivity should be finite, smooth, 
and independent of t, leading to 



lim T<±(x) — sqx 



= »„*(l-«.y)/2 



(33) 



where sq is a universal number. Thus, at T = T c , 

an = {SkS^B/HuY 1 -^' 2 . (34) 

A similar limit can also be taken for the conjugate fields; 
we show only the results for H ex : 



CT fl = (s k s Q ){H k n) 



%' 2 



Hq + H e 



(35) 



The asymptotic behavior of £_(x) as T — * T m (B > 0) 
is of particular interest when the transition is continuous. 
The multicritical description 133] involves a crossover 
from XY to glass fluctuations |5[-]^jL9| . Although glass 
fluctuations dominate near T m {B > 0), the XY scaling 
fomula d32|), is still appropriate. Approaching the tran- 
sition, we find 



lim oc 



\x £ m ) 



X *^-"Ti 

X Xfi 



(36) 



where x m is the same universal constant as in Section IV. 
The exponent u> is related to glass, not XY fluctuations. 
In Ref. 1 19 1 it has been argued that id — v g {z g — I), in 
analogy with Eq. (pi]) ) where v g and z g are glass expo- 
nents. 



IX. CONCLUSIONS 

In this paper, we have presented the basic scaling the- 
ory of the intermediate 3D XY transition. The theory 
is quite general; it involves only the assumption that 
B ~ (Length) , which is deduced from minimal cou- 
pling. We stress that the theory should apply to all 
strongly type-II superconductors, including both high-T c 
and low-T c varieties. In particular, we have considered 
the case that the finite field transition T m (B) is contin- 
uous, although a similar analysis can be applied in the 
case of a first order melting transition p2[ . 



The scaling results can be summarized by noting that 
the H-T or H cx -T superconducting phase diagrams in- 
volve two regimes, as shown in Fig. 1. In the outer region, 
scaling is the same for each of the fields B, H, and H cx , 
up to very small correction terms. To obtain results for 
the different fields, we consider the equations involving 
only B [for example, @, (|), @, @, and ©], then 
apply the following substitutions: 



B 



H 



He 



(37) 



In the inner region, scaling behaviors differ for B, H, 
and H cx . This is a nontrivial consequence of the conju- 
gate nature of B and H, in the thermodynamic sense. 
Observation of the inner region should therefore be re- 
garded as a more stringent test of 3D XY scaling. How- 
ever, care must be taken to distinguish inner scaling be- 
havior from inverted XY behavior. The dimensions of 
the inner region, in the H-T plane, are given by the char- 
acteristic field and temperature scales, Hq and t\T c . In 
the H ex -T plane, these become Hq and t\T c . 

The high-T c oxide materials are natural candidates for 
observing both the inner and outer 3D XY scaling be- 
haviors, due to their strongly type-II character, and the 
prevalence of vortex fluctuations near the superconduct- 
ing transition. However, recent experiments demonstrate 
behavior consistent only with the outer region. To under- 
stand this, it is helpful to obtain estimates for Hq, H e g, 
and t\. We can make use of two experimental analyses. 
The first (I), by Hubbard et al. [ ^2| , involves a series of 
YBa2Cu307_<5 single crystals with varying 8. The sec- 
ond (II), by Moloni et al. [Q, uses a similar series of thin 
films. Both cases involve samples with T c ~ 77 K. Us- 
ing the appropriate demagnetizing factors (n ~ 0.84 for 
sample I, and n ~ 0.992 for sample II), we can make our 
estimates. We will refer to the geometry-independent 
properties of each sample (without tildes) as "intrinsic." 
These intrinsic quantities are assumed to be the same for 
both of the 77 K samples. 

In Hubbard et al., the fluctuation magnetization was 
plotted for T ~ T c , as a function of the field H cx . We 
can then use Eq. ( p7j ) to find Hq, with the assumption 
Cl c ~ 1. This gives Ho ~ 5 x 10~ 7 Oe for sample I, corre- 
sponding to the intrinsic result Ho — 2 x 10~ 5 Oe. The 
estimate for sample II then becomes Ho — 1 x 10~ 9 Oe. 

The relative temperature scale t\ can be estimated by 
assuming that along the superconducting transition line, 
H cx ~ Hq when |f| = t\. From Eq. @, this gives 
l x ~ {H /2H*) 1 / 2u *y, where H* = Q m x m H k is another 
characteristic field of experimental significance. For sam- 
ple II, it was found that H* ~ 19 T. We can then estimate 
t\ ~ 1 x 10~ 9 for sample I and t\ ~ lx 10~ n for sample II, 
with the corresponding intrinsic result, t\ ~ 2 x 10~ 8 . 
Note that the width of the inner scaling region of sam- 
ple I becomes t{T c ~ 9 x 10~ 8 K. 
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The estimates given above for Hq and t\ pertain to 
underdoped high-T c cuprates. For the case of opti- 
mally doped cuprates, estimates for Hq and t\ should 
be somewhat smaller. Estimates for low-T c superconduc- 
tors should be even smaller. In each case, the smallness 
of the inner scaling region reflects the weak diagmag- 
netic response of fluctuations above T m . Thus the inner 
region is probably experimentally inaccessible in many 
cases, due to sample inhomogeneities. However, strong 
anisotropy effects may improve the situation. As dis- 
cussed in Sec. II, an assumed relation (fk/jT c ) 2 oc i?| 
leads to H oc 7 2 T ( ?, showing that H (and ii) may be 
greatly enhanced in very anisotropic samples. 

It is also possible to estimate of the temperature width, 
tj nv , of the inverted 3D XY scaling region. This may 
be compared to the size of the inner scaling region, t\. 
Although neither the intermediate nor inverted 3D XY 
models has been solved exactly, the crossover may be es- 
timated using the intermediate scaling theory: it is the 
point where our assumption of a uniform magnetic field 
breaks down. (See Sec. II.) Approaching the transition, 
the zero-field screening length A = \ \t\~"*y /2 , 

diverges more slowly than the 3D XY critical correla- 
tion length, £ = £oKI -I/x!/ - We expect a crossover to 
occur when k — A/£ ~ l/\/2. The length scale £o of 
3D XY fluctuations is not well known experimentally; 
however, a rough estimate can be obtained by using the 
bare Ginzburg-Landau parameter, which is easily ob- 
tained for samples such as I and II: k — A /£o — 100- 
250. We finally obtain t inv ~ 2-40 x 10" 8 . In this rough 
estimate, the inner and inverted scaling regions are of 
comparable size, and must both be considered in order 
to correctly interpret fluctuations effects very near the 
3D XY critical point. 

In contrast to the elusive inner scaling region, the outer 
region, or behavior consistent with it, is readily observed 
over a wide temperature range: \t\ < 0.5 0, which is on 
the order of 10 9 -10 11 times t\. The difference between the 
two XY temperature scales is striking. However, we note 
that there is no reason why they should be related. It 
remains an outstanding theoretical problem to provide 
estimates for the relevant temperature and field scales, 
from microscopic considerations. 
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